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Abstract. A didactic description of charge and spin equilibrium currents on 
mesoscopic rings in the presence of Spin-Orbit interaction is presented. Emphasis 
is made on the non trivial construction of the correct Hamiltonian in polar coordi- 
nates, the calculation of eigenvalues and eigenfunctions and the symmetries of the 
ground state properties. Spin currents are derived following an intuitive defini- 
O , tion and then a more thorough derivation is built upon the canonical Lagrangian 

O i formulation that emphasizes the SU(2) gauge structure of the transport problem 

of spin 1/2 fermions in spin-orbit active media. The quantization conditions that 
follow from the constraint of single-valued Pauli spinors are also discussed. The 
CN ' targeted students are those of a graduate Condensed Matter Physics course. 

QiQ Resume. On donne une description didactique des courants de charge et de 

1^^, | spin dans un anneau mesoscopique en presence de couplage spin-orbite. L'accent 

est mis sur la construction de l'hamiltonien, le calcul des valeurs propres et 
des vecteurs propres et les proprietes liees a la symetrie par renversement tem- 
porel. Le courant de spin est calcule a partir d'une definition intuitive, puis 
l/~) ■ une demonstration plus approfondie est construite a partir de la formulation la- 

f— *^ ■ grangienne qui met en evidence la structure de jauge SU(2). Les conditions de 

^"■s ' quantification qui decoulent des contraintes de raccordement des spineurs de Pauli 

sont egalement discutees. Le public vise est celui d'un cours avance de physique 
de la maticre condensee. 
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1. Introduction 

The study of spin transport properties in materials, especially in semi-conductors 
where two-dimensional electron gases (2DEG) can be fabricated, has become an 
important field of research both theoretically and experimentally in the last decade. 
One of the key concept there is the Spin-Orbit (SO) interaction which allows for spin 
manipulations. Mesoscopic rings are very simple laboratories for the investigation of 
various quantum effects PQ in the presence of spin-orbit interactions. They are the 
subject of intensive interest, both for fundamental reasons and potential applications 
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in spintronics devices. Let us only mention a few fundamental studies in connection 
with phase effects, such as Aharonov-Casher effect [2]l3] or Sagnac phase shifts [4], spin 
interferences and spin filtering |5-9] or in the fashionable subject of graphene band 
structure [10]. On the more applied side, we can cite spin manipulations [TOfM] . 
studies of spin related transport properties [TSHTT] and persistent currents [TM2TJ] . 
spin filtering in arrays of SO quantum rings |21j that operate on principles of spin 
interferometry analogous to their optical counterparts [22j and concepts for the 
elaboration of qubits gates ;23 a . 

Many of the concepts involved in spin transport are rooted in very basic quantum 
mechanics but in contexts that may be unfamiliar for the beginner, such as multiply 
connected geometries like closed loops or rings. This paper will serve to sort out such 
uncommon applications opening the scope of understanding forefront topics such as 
spintronics. 

In this paper, we present an analysis of spin transport in a mesoscopic ring. 
The exposition is suitable for a course in condensed matter physics at the graduate 
level. We carefully present the details of the non trivial construction of the correct 
Hamiltonian in polar coordinates, the calculation of eigenvalues and eigenfunctions 
and the symmetries of the ground state properties with respect to time reversal 
symmetry. Charge and spin currents are derived following an intuitive and heuristic 
definition and then a more thorough derivation is built upon the canonical Lagrangian 
formulation that emphasizes the SU(2) gauge structure of the transport problem. The 
path integral approach and the quantization conditions that follow from the constraint 
of single-valued Pauli spinors are also discussed for the benefit of a full picture that 
will be valuable to the student. 

The Spin-Orbit interaction finds its origin in the Pauli equation, which follows 
from the non-relativistic limit of the Dirac equation, 
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Here we consider electrons, e = — |e|. Ordinary vectors are denoted by arrows and bold 
face characters are used for 2x2 matrices. The first term in the first line corresponds to 
the usual Schrodinger equation including the kinetic energy with a mininal coupling to 
the electromagnetic gauge field A and a scalar potential contribution — e<f>. The second 
term in the first line describes the first relativistic correction to the kinetic energy and 
the Darwin term, where E is the electric field and c the speed of light. These first two 
terms are proportional to the 2x2 identity matrix in spin space 12x2- The second line 
comprises explicitly spin-dependent terms, first the Zeeman interaction [24] where B 
is the magnetic field and er is the vector of Pauli matrices [21] and the second term 
is the Spin-Orbit interaction, written with the minimal coupling to the gauge vector. 
We have assumed a static potential so that the rotor of the electric field is absent and 
the Spin-Orbit interaction is limited to the term mentioned here |25j . 

Forgetting about the magnetic field, the Spin- Orbit interaction has a simple 
interpretation in terms of the interaction of the spin magnetic moment of the particles, 
(here supposed to be electrons of Lande factor [26] g — 2 and spin s = \hcr), 

P 1 = 92m^ = — 2m ^ with the magnetic field produced by all external moving 
charges in the electron rest frame. Assuming a uniform electric field E acting on the 
moving electrons, the SO contribution to the Hamiltonian results from the interaction 
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— /I-B regt £ r of p, with the effective magnetic field experienced by the particles in their 
rest frame. In the case of a Lorentz change of reference frame, one has 

B r est fr. = Jr(-#) X E = -(mc 2 )" 1 ^ X E), (2) 

where v and p refer to the dynamical variables of the electron. This expression is 
corrected by a similar contribution with a factor — \ due to the Thomas precession in 
the case of a closed orbit. Such a term appears when the proper Lorentz transformation 
for the fields is considered [26) . In such a way, the resulting SO interaction is usually 
written (e.g. in atoms) 

tfso=-^-(pxE), 

_ N lWr) rl , (3) 



m 2 c 2 r dr 

were we have substituted a spherically symmetric potential, E = — Vcf> — - ^ r and 
the interaction is proportional to s • L. It is clear now why, in the context of atoms, 
the interaction is christened the Spin-Orbit interaction i.e. L pertains to the angular 
momentum of the orbit while s pertains to intrinsic angular momentum, the spin. 

The prefactor in this expression depends on the details of the problem (Lande 
factor, effective mass, . . . ) its sign even depends on the nature of the particles 
involved. In semiconductor physics, one usually introduces the notations a or j3 for this 
coefficient, which is determined perturbatively in the k -p theory by a matrix element 
in the Bloch wave functions basis. In the case of a two-dimensional electron gas with a 
gate voltage applied perpendicular to the 2d sample (with a non-symmetric confining 
potential generating a space inversion assymmetry (SIA)), this term is known as the 
Rashba SO interaction l27ll28l. 



V Rashba = const. 9 ■ (E X p) = a(er y p x - cr x p y ). (4) 

Note that the Rashba SO amplitude can be tuned experimentally using a gate voltage, 
since the prefactor a is proportional to the electric field (see figure [1} . 

Gate Voltage 



E 



2D Electron gas 



Figure 1. The typical Rashba situation with an electric field, created by a gate 
voltage, perpendicular to a ring fabricated on a two dimensional electron gas. 

The calculation of matrix elements leading to a is generally a hard task and one 
usually uses phenomenological expressions compatible with the crystal symmetries or 
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computes them directly from experiment |29j . The Rashba interaction applies in the 
case of SIA. When there is bulk inversion asymmetry (BIA) i.e. the crystal unit cell 
lacks inversion symmetry, we have the Dresselhaus 'flavor' of the SO interaction in 3d 
systems [30] . 

Vr> j3d = const. k x (k 2 - k 2 z )cr x + c.p. (5) 

where c.p. stands for cyclic permutations. In the case of electrons confined in two 
dimensions, the expectation value along the confinement, quantized dimension should 
be considered. If that direction is z, then (k z ) ~ 0, {k 2 ) ~ (n/£) 2 , I being the typical 
confinement length or spatial width of the confinement potential. The Dresselhaus SO 
interaction thus takes the simple form, here written in a notation closer to equation l|4|). 

V Dresselhaus = f3{&xPx ~ CyPy), (6) 

where we neglect cubic terms in k and the average values computed are lumped 
in the definition of (5. For more details on SO interactions in semi-conductors, see 
Refs. [30H32]. 



2. Derivation of the ring Hamiltonian 

2.1. The argument of Meijer, Morpurgo and Klapwijk for the Rashba SO interaction 

To make the discussion simpler, we neglect any magnetic field effect, which is well 
discussed in the literature. Here we discuss the non-trivial point of designing a 
Hcrmitian Hamiltonian when it is not in cartesian coordinates. We explain in detail 
such construction so the student can understand the pitfalls generally avoided in 
quantum mechanics courses because of a preferential system of coordinates. Many 
researchers overlooked this hermiticity problem until the recent paper by Meijer, 
Morpurgo and Klapwijk [33] - 

The classical argument in reference (33] is essentially the following. We consider a 
2DEG, neglecting interactions between electrons and add the Rashba SO interaction 
in cylindrical coordinates (p, ip, z), 

H 2d = -^(d 2 p+p -% + p - 2 dl)l 2x2 

+ap~ 1 (cr x costp + (T y siinp)ihd v + a(cr x shop — cr y cos <p)ihd p . 

The Rashba SO interaction in the second line comes from equation Q with the 
definitions p x = —ihd x and p y = —ihd y after the substitutions d x = — (p)^ 1 simpdp + 
cos ip dp and d y — {p)~ l cos ip d v + smtp d p . Fixing the radial distance p to the ring 
radius a and neglecting the radial derivatives leads to a "Id" Hamiltonian (called H'i^ 
for the moment) for the ring, 

fr 2 

H'id = ~{id v ) 2 + aha (<r x cos<p + er y shx(p)id v . (8) 

In [33] . the authors note that the last term in (|5]) is non-Hermitian, since 

(F\(<r x coB<p + tr y Bm<p)iB v \G)' = (J^dfW K) ( ie £ d(p Ze "^) (gj 

= (G\(cr x cosip + a-y sin(p)id v \F) — i(G\(cr x siny> — cr v cosy)|F), (9) 
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where \F) and \G) are Pauli spinors, and the integrated terms vanish if we impose 
single-valued spinors, i.e. (ip + 27r\F) = (<p\F). From the property given in equation 
(|9|). we may infer that a Hermitian operator could be formed by the combination 
V = aTia W where 

W = {<T X cosy + <r y sa\Lp)id v — iA(cr x simp — cr y cosy), (10) 

(11) 



since 



(F\W\G)* = (G\ (cr x cos p + (T y sin tp)id v \F) 

-j(l - A*)(G\(cr x siny - er y cosy)|F), 



and 



(12) 



(G\W\F) = (G\ (a- x cos <p + (Ty sin p)id v \F) 
—iA{G\{cr x siny — cr y cos (p)\F). 

The r.h.s. of equations (flT]) and (fT2"j) are equal provided that A* = A = i. This value 
in equation (|10[) gives the correct expression that renders the Rashba SO interaction 
on the ring Hermitian. 

The method used in Ref. 33 to build hcrmiticity into the SO Hamiltonian, 
consists in the introduction of a confining potential V(p) in equation ([7]) in order 
to localize the particle on a circle of radius a. If this potential is steep enough, the 
particles will lie in the ground state i?o (p) of 

H p = KE + V(p) = -A- (d 2 + p-'dp) + V(p), (13) 

where KE is the contribution of the radial motion to the kinetic energy and the correct 
Id Hamiltonian (H^f /lf)a ) is defined by the projection, on this ground state, of the 
angular dependence of H^d, 

H r i, = <^o|H 2d - KEl2x2l.R0), (14) 

where \Ro) is extended to a two-component object. Note the very peculiar situation 
in polar coordinates, KE does not simply coincide with p' \/2m (see e.g. Ref. [33] for 
a discussion on this unrecognized point). The calculation is now made specific by 
considering a circular Harmonic potential V(p) = ^rnu 2 (p — a) 2 and by the explicit 

,-1 n „A II3-\P) I 



calculation of the matrix elements (Ra\p 1 \Ro) — a 1 and (Ro\d p \Ro) = —(2a) 1 . 
The resulting Hamiltonian H„ Q ^ /if)a then follows from the substitutions p — » a, 



Rashba 

d p -> -(2a)- 1 in H 2d + £ (d 2 p + p- l d p ). Thus 

rin ft 2 

H « h = 2W ( ^ )2l2x2 (15) 

+aha~ 1 (cr x cos(p + <r y sin •p)id Lp — iah(2a)~ 1 (cr x shxcp — cr y cosy). 

It can be checked that the Hamiltonian is now Hermitian, since it coincides with the 
choice A* = A = |, according to the procedure outlined in Eqs. (|10II12[) . 

2.2. Symmetrization of the original Hamiltonian 

Although the previous derivation is of course correct, it is instructive to propose 
another, probably more straightforward derivation. The SO interaction in cylindrical 
coordinates with a perpendicular electric field, a = ae z , takes the form 

<x • (a xp) = -aio-pp^ -o-^pp), (16) 

with (T p = cr x cos ip + cr y sin ip and cr v = — <r x sin ip + cr y cos tp treated for the moment 
as classical variables (i.e. we forget about possible commutation problems). In order 
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to make the corresponding expression Hermitian in Quantum Mechanics, one has 
to use the correspondence principle, i.e. add inequivalent orders of non-commuting 
operators, and symmetrize the "classical" expression, using also p p = —ih(d p — (2p)^ 1 ) 
and p v = —iha~ 1 d ip . Let us also note that when the electrons are confined on a ring, 
then (p p ) = (R \p p \Ro) — (and not (Ro\d p \Ro) — as it was done to get equation ([5])) 
and p = a. This leads to 

Y R^hba = -2 a {°>'M = i^a- 1 ((T p d v + -d v <r p ), (17) 

where {<T p ,p v } denotes the anticommutator cr p p v + p v cr p . Using d v cr p — cr v , 
equation (fT5j) eventually follows, 

rin S 2 1 

H SaL,a = ^^^) 2 l 2 x2 + lhaa \(T p d v + -0>). (18) 

Let us emphasize on the fact that the condition (p p ) = is equivalent to Meijer et al's 
derivation, since it implies —ih(Ro\d p — (2p)~ 1 \Ra) = with a radial wave function 
strongly localized on p ~ a. Our proposed derivation, nevertheless, has the merit of 
being simpler and more general, since it does not rely on any particular assumption 
on the confining potential (except the fact that it enforces the particles to move along 
a ring). 

2.3. Hamiltonian in the rotated spin basis 

In the previous sections, the Rashba Hamiltonian (|18|) . initially expressed in Cartesian 
coordinates, has been rewritten in polar coordinates. To preserve the hermiticity of the 
Hamiltonian, it has been shown that an extra term should be added if the expression 
is not properly symmetrized prior to quantization. In the following, we give a more 
physical picture that gives rise to this extra term. 

We first start by noting that the change of coordinates (x,y) — > (p,<p) induces 
a rotation by an angle ip of the local frame: (e Xl e y ) — > (e*p,e v ). All vector fields, 
including a are now expressed in this local frame. However, no such rotation has been 
applied to the spinor basis which thus remains the same at any point of the circle. 
The notations a p and cr v may thus be misleading. Let us introduce such a rotation 
of the spinor basis: 

|F')=e 4 f-|F) = (cos|l 2x2 +zsin|^)|F). (19) 

Since the rotation axis is (Oz), the quantification axis is not changed during the 
transformation. The two spin components simply get a phase 

|f> = e**|t>, |V) = e-**\l). (20) 

In the new local spinor basis, the Hamiltonian becomes 

H^e'^He-'^. (21) 

Using the identity cr p er z = —icr v , the correct Hermitian Rashba Hamiltonian can be 
cast as 

Htfhia = ^ 2 {^?^+^a-^ p (d v + i ff ,). (22) 

Under the change of spinor basis, the new rotated Hamiltonian becomes 

e'S-'H^e-**" = ^ (id v Ux2 + \rr z f + l haa- 1 <r x d^ (23) 

since cr' = e l ~ a ~*cr p e~ l ~' Tz — cr x . Interestingly, after rotation of the spinor basis the 
potential energy recovers a form ex cr'p v similar to the original Rashba Hamiltonian 
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before symmetrization in (j 16[) . From the physical point of view, one indeed expects 
the potential energy to be the same in all local orthonormal frames. This implies 
that, starting from the original Rashba Hamiltonian, the change of coordinates should 
be accompanied by a change of the local vector and spinor frames. The spin-orbit 
term then keeps the same form. When expressed in the original spin basis, the 
additional term needed for symmetrization appears naturally. In contradistinction 
to the potential energy, additional terms may appear in the kinetic energy when one 
changes the spinor basis. They are interpreted as inertial forces due to the fact that the 
local frame is not inertial. This is the origin of the extra term \<* z m the Hamiltonian 
(|23[) (note the appearance of the z component of the total angular momentum, L z +s,z). 
When the particle moves along the circle, the local spinor basis turns by the same 
angle. Consider for example a constant spinor. To accommodate the change of local 
basis along the circle, it should have (^-dependent spin-components. Since no kinetic 
energy is associated to the spin, a term involving derivatives d v would produce a term 
which is precisely cancelled by the additional contribution \<r z - 

3. Eigenenergies and Eigenvectors on a ring 

3.1. The Rashba SO interaction 

Once the correct Hermitian Hamiltonian has been written in polar coordinates we can 
confidently derive the eigenvalues and eigenvectors for the strictly one dimensional 
ring. This section will benefit the student by explicitly obtaining the eigenfunctions 
and eigenvalues of the system, from which we can obtain all observables or actually 
measurable quantities in experiments. We will also make some basic symmetry 
considerations that will illuminate details of why the wavefunction and energies take 
their particular form. 

First we rewrite the Hamiltonian in a clever way in order to arrive rapidly at the 
eigenvalues and introduce physical insight by pointing out the existence of a gauge 
field associated to topological phenomena connected to spin transport. An alternative 
form of the Hamiltonian (fT8|) was given in Ref . [35] . We observe that 

/ maa \ 2 . , 2 2maa maa [maa\ 2 

U<9 v l2x2H -r— a-pj = [id v ) l 2x2 H t — <r p {io ip ) + — — ier v +y— — j l 2x2 .(24) 

The factor of two in the second term on the right accounts for the appropriate action 
of an derivative operator. Such form allows one to recast the Hamiltonian as 

/ „ maa \ 2 /maa\ 2 

(«c^l 2x2 + —— <7 p J - {—) l 2 x 2 ■ (25) 

In order to find the eigenstates of the Hamiltonian (|2"5j) . one first solves the 

eigenvalue equation 

/ „ maa \ 

[id v l 2 x2 + —f-vp) * = £ *> (26) 
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It follows that the general form for the cigenfuctions should be 



(27) 



/4, e -^/ 2 \ 

*t-e^W^ e elW2 ), (28) 
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Figure 2. Energies of electrons on a ring with spin-orbit interaction as a function 
of the SO interaction strength a. The energy is in units of h 2 /ma? while a is in 
units of h/ma. Indicated are the eigenvalue labels for the first five levels (note 
the ground state is fourfold degenerate). The graph shows how the limit a = 
is fourfold degenerate as expected since the time inversion symmetry turns into 
independent space and spin inversion symmetry. Arrows indicate change in lowest 
energy level. 



where n is the main quantum number and we account for the two eigenvalues that 
will arise, s = ±, for each wave propagation direction A = ±. So, assuming such a 
wavefunction we obtain the corresponding eigenvalues 



8 I /maa\ 
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(29) 



2 V V h 

where n is a half odd integer as will be seen below when we address the wavefunctions. 
The eigenvalues of (|23|) then simply follow 

2 



E n.s - 



2ma 2 



As , 

Y^ 1 



V 1+4 hr) 



(30) 



Note that left and right propagating waves with the same s-index are not degenerate 
but time reversal symmetry is satisfied i.e. E* s = E~ _ s (simultaneous change of A 
and s). This symmetry reflects the fact that the SO interaction is not space inversion 
symmetric, but only time reversal symmetric. When the spin-orbit interaction is 
absent (a = 0), space and spin inversion symmetry is recovered and energies of both 
clockwise and counterclockwise modes are degenerate independent of the spin label. 
Figure [H shows the energy levels as a function of the SO strength. The ordering of 
the levels are indicated according to the spin orientation and sense of the current. 
The free electron on a ring is recovered for all levels as one can see from the fourfold 
degeneracy at a = 0, rendering all possible combinations of the values of A, s for a 
fixed integer value of TV in the relation E^ s = h 2 (n — \s/2) 2 /2ma 2 — h 2 N 2 /2m,a 2 
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Table 1. Energies (in units of H 2 /2ma 2 ) for the first three levels and their 
degeneracies according to values of the quantum numbers n, A and s in the limit of 
zero spin orbit, and the corresponding values of the free electron integer quantum 
number N. Also shown is how degeneracies occur when the SO interaction is 
turned on leaving only time reversal symmetry. 
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(with n half odd integers) . Table 1 shows the energies for the first three levels in the 
limit of a = and how degeneracies are broken when the SO interaction is turned on. 
Now we compute the eigenfunctions. In order to deal with the subtleties 
associated with the change in sign of spin and propagation direction we explicitly 
do the A = + case. Using equations (f2~T|) and (|2"5|) we find from the secular equation 

„— (^-lW, s , (31) 

2maa \cos9 I 



B+.s — 



where cos (9 = l/\/l + A(maa/h) 2 . In order to conform to a canonical spinor we 
choose A+.+ — cos 6/2, and thus B+,+ = sin 9/2 for normalization, so we have 

h /1-COS0A 9 9 

cos - = sin -, (32) 



2maa \ cos 9 J 2 2' 

from which we have the condition 

2maa 
tan0=— — . 33 

h 

The choice for the second eigenfunction is A+ t _ = —sin 9/2, and doing the same 
exercise we arrive at the two eigenspinors 



Te"^ /2 \ T + iw f- sin le-^/ 2 



xh+ _ p my[ UJt, 2 e \ fl>+ - pimp I bm 2 e | fiA\ 



The corresponding eigenfunctions for A = — are 

/ cos ^e"^/ 2 \ ( - sin ^e - ^/ 2 A 

^ + = e ^(s?nfeW 2 J' ^- = e ^( c^fe^ ) ^ 

Using the time reversal operator for spin 1/2 particles = —icr y K where <r y is the 
corresponding Pauli matrix and K is the conjugation operator [36] . one can readily 
show that 1 &n+ = ®^n- an d ^n- = ®*n+- ^ s one would expect from the 



Icharge = - 2_^ f,^ (36) 
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time reversal invariance ^~ + is degenerate with \E r „ _ and \t^_ with \E r „ , . With 
these properties in mind the student can guess that expectation values taken with 
corresponding time reversed wavefunctions must be added together since there energies 
are degenerate. This fact will have surprising consequences in the next section. 

4. Charge and spin currents in the ground state 

4-.1. Direct calculation 

Having found the eigenfunctions, one can now compute equilibrium properties such 
as persistent charge and spin currents on the ring. The student might be accustomed 
to expecting a current only when an external force is applied to the system, such as 
a potential difference. Nevertheless we will show that when symmetry breaking fields 
exist that can do no work, such as a magnetic field, currents can exist as equilibrium 
properties borne from the nature of the wavefunctions. In a sense the system can 
distinguish between clockwise moving current and a counterclockwise moving current, 
one of them corresponding to a lower energy. These currents are long lived in the 
absence of perturbations. This is very counterintuitive for the classical line of thought. 
There are charge currents in equilibrium only when time reversal symmetry 
breaking perturbations are present, such as a magnetic field [T1137] . for an experimental 
review see [38] and references therein. The charge currents can be derived directly 
from the ground state energy as 

where $ is the magnetic flux and the A quantum number is chosen by the direction 
of the magnetic field. An alternative form of computing the currents is by using the 
definition of the current operator. These current densities are fields usually defined 
from a continuity equation rather than from matrix elements of some operators. In 
order to generalize the classical definition of the current density j(r) = n(r)ev, we 
take the expectation value of the velocity operator 

Jcharge = * f eV *. (37) 

where e is the electron charge and v is the velocity operator. Note that the current is 
defined in such a way that its dimension is a charge times a velocity. In the presence 
of the spin-orbit interaction, the velocity operator is not simply p/m. There arises 
what is known as an additional anomalous velocity term. We start from the quantum 
mechanical definition of the velocity (which is now a 2 x 2 matrix due to the presence 
of the spin-dependent terms in the Hamiltonian) v = |[H,r|. Since we are interested 
in the azimuthal velocity component ip, it is more convenient to calculate 

v ^ = -r [H, <p] = - <9 v l-2x2 -acT p . (38) 

h ima 

rather than working in Cartesian coordinates. Note that the Hamiltonian takes 

a simple form when expressed in terms of the velocity: H = -^rniy 1 — a 2 l_2x2)- 

Here, we do not consider symmetry breaking magnetic fields, so there can be no 

persistent charge currents on the ring. If we compute the charge currents for the 

fourfold degenerate lowest energy levels from which one forms a totally symmetric 

linear combination, one gets 

" |(*: 1/a ,+) t v».*: 1 / ai+ - |(^ 1/2 ,-) f v^+ 1/2j _ (39) 
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(40) 



provided the SO strength is within the arrows in Fig. [5J since outside the interval the 
computation involves a different two-fold degenerate level. The result above applies 
to any normalized linear combination of the fourfold ground states satisfying the 
time-reversal symmetry. Charge current cancellation happens level by level and for 
all a values, rendering them zero as expected. The effect of an external magnetic 
field (playing the role of a time-reversal symmetry breaking field) and the ensuing 
charge persistent currents in the presence of a simple scalar potential and the spin- 
orbit interaction were examined in ref. [18] . The small scalar potential breaks the 
degeneracy at the edge of the Brillouin zone that limits the charge persistent currents 
to a maximum value for each band. 

On the other hand, there can exist currents that don't break time reversal 
symmetry i.e. spin currents. In order to make an explicit calculation of the current, 
here we sidestep the consideration of a scalar potential and use the empty lattice 
approximation [39] . The spin current density, denoted as J a in order not to get 
confused with the charge current density J, follows from the same approach as above. 
Note that the spin current J a is a tensor (two indices) while J is a vector. We define 
now a local "spin velocity" operator, properly symmetrized 



J a = ^{v, S a }V, 



(41) 



where s a = hcr a /2. The reader should note that the time reversal operation applied 
to such an operator (in the absence of a magnetic field) reverses both the velocity and 
the spin, so the current is not changed by such an operation, so this intuitive definition 
materializes our earlier expectation. 

The spin current along the ring will again involve the previously derived azimuthal 



velocity operator 
current is given by 



The general contribution for any single state involved in the 
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As I 



(42) 



where a — x,y, z and we have used notation from Eq. (|28[) . Adding contributions from 
the four degenerate lowest lying levels (again the totally symmetric linear combination) 
with the spin orientation in the z direction, one arrives at 
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(44) 
(45) 
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Figure 3. Spin current in units of h 2 /4ma from the lowest lying degenerate level 
indicated in Fig[2]as a function of a. Between the arrows in Fig[2]the lowest level 
is fourfold degenerate so there are four terms in the spin current. Outside that 
range, other two-fold degenerate level take over causing a jump in the spin current 
as indicated. The graph shows how the limit a = indicates zero spin current 
because up down spin symmetry along with inversion symmetry are restored. 



Such an expression only accounts for a values in the range between the arrows in 
FiglH Outside this range the energies E^, % , E~, 2 _ are the lowest energies, so 
we must compute the expectation of the spin current with such two-fold degenerate 
eigenfunctions. The full spin current for the lowest lying levels is depicted in Fig [3] 
One can see from the figure, that at a = the spin currents vanish, as required by the 
recovery of the inversion symmetry (fc — > —k) and, independently, the spin inversion 
symmetry that were intertwined in the presence of the SO interaction. In the presence 
of a small scalar potential on the ring, the level crossing in FigJ5] gives way to level 
repulsion. The abrupt transition in dashed lines in Fig|3] will then be gradual and 
rounded |37j . Persistent spin currents have yet to measured due to the lack of an 
appropriate probe. 

The spin currents J* ,v can also be readily computed along the lines above. In 
contrast to the spin current J* , such currents for each of the degenerate states of a 
particular level, are not real valued. Nevertheless, the summation of all the degenerate 
contribution renders a real spin currents which for the ground state is (in the spin-orbit 
strength range shown in Figj2j) 



$v = -rfra(cos6>- l)cos<p, 
jv = -ha(cos9 — l)sin<^. 



(46) 



Such expressions have a cp dependence due to the precession of the spin around the 
z axis. This means that a particular polarization rotates going through zero and 
then changing sign, making the whole spin current change sign while the wavevector 
is constant. Note that the spin persistent currents vanish correctly for a = as 
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expected. Figure [4] shows the angular dependence for both spin current components 
for a value of a within the range shown in Fig 12] 



0.3 r 




Figure 4. Spin currents for the x,y polarizations in units of fi 2 /4ma from 
the lowest lying degenerate levels indicated in Fig[2] as a function of tp and for 
a = „ , Between the arrows in Fig[2] the lowest level is fourfold degenerate so 
there are four terms in each of the spin currents. 

As with charge persistent currents one has to add all occupied levels. Observing 
figure [5] one sees that the first two levels add constructively while the third level 
subtracts from the previous two levels. One must then carefully add the corresponding 
currents taking into account level crossings. 



4-2. The Dresselhaus SO interaction 

The Dresselhaus Hamiltonian of equation (J6j> is another form of SO interaction. 
Using either of the procedures shown in section 2, keeping only p v terms, we express 
CxPx — cfyPy as — <r x sin tppp — cr y cos pp v that we then symmetrize to get 



yrmg 

Dresselhaus 



■\(3(a x {sinp,p v } + a v {cosp,p v }) 



(47) 



= ih[3a 1 [(cr x simp + a y cos tp)d v + ^(cr x cos p — cr y simp)}, 
or we define the combination 

W = (<r x simp + (T y cos ip)id v + iA(a x cos tp — cr y simp), (48) 

and form the Hermitian conjugate 

W T = (<t x simp + cr y cosip)id v + i(l — A*)(cr x cosip — <J y wmp), (49) 

where hermiticity requires A* = A = i 

The method of Ref. [33] applied to the Dresselhaus SO interaction would require 
to express d x and d y in equation ([6]) in terms of d p and d v , then to substitute p by a 
and dp by —(2a) -1 . It obviously leads to the same expression. 
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Now, completing the square as we have done in the Rashba case, we get for the 
Dresselhaus SO interaction the more compact expression 



•erring _ "- 

Dresselhaus Ornn^ 



io v l2x2-\ r— (cr x sin <p + cr y cos <p) \ - — r— 1-2x2 



.(50) 



EL 



2ma 



This expression is very similar to the Rashba Hamiltonian and only slight modifications 
are needed to calculate the spin current along the ring. The energy levels are now 
given by 

-^y/l + A(ma(3/h)A - (ma(3/h) 2 ] , (51) 

i.e. they are obtained from the Rashba eigenenergies via the substitution a — > ft and 
time reversal symmetry is preserved. The corresponding eigenspinors become 

<+ = etnv ( -uM^ ) ' *-- = e " lv ( -i ™l^ ) (52) 

sin|e^/ 2 \ .„._ _,-„,„/ -cos |e^/ 2 



*«.+ = e ~ mv -Zk.-^ .**- - ^ SZ L-^ ) ( 53 ) 



—i sm ttC 



where now cos 9 = 1/yl + 4:(maf3 / h) 2 . When the SO coupling (3 is not too strong, 
we recover the structure of a four-fold degenarate ground state and the spin current 
polarized in the z direction takes the form 

^ = ^(1 -cos.). (54) 

Notice that the computed z component of the spin current has the opposite sign to 
that of the Rashba contribution (see (|45p). so that given equal SO strengths for both 
type of interactions the two effects compete to produce a vanishing net spin current. 

4-3. Spin current and the non-Abelian gauge formalism 

The computed spin currents in the previous sections where based on intuitive 
extensions of the very familiar charge current definition, i.e. the physical quantity 
being transported times the velocity, along with the necessary symmctrization when 
the definition involves two non-commuting operators. In this section we make the 
student realize that a formal, well grounded approach can be used to define the currents 
based on regarding them as conserved quantities following Noether's theorem. We will 
depart from a general Lagrangian that will be identified with the corresponding ring 
Hamiltonian, derived above. New tensor fields will arise, analogous to the gauge 
potential in electromagnetism but of a non-abelian nature, that will by canonical 
relations, give us the full expressions of the currents in the theory. 

The spin current calculated above is the equivalent of the paramagnetic current 
in the case of U{1) gauge theory (see e.g. [5]). One can ask about a possible equivalent 
to the diamagnetic contribution also (such a contribution is called color diamagnetism 
by Tokatly [40,41 V The Lagrangian approach is very convenient to investigate this 
question. To the Hamiltonian given in equation (|25|) . we associate a Lagrangian density 
£ defined according to 



L — I adip C, 
= <*|ififtl axa -H|*), 
= [ ad(p[ih&& -tfHV]. (55) 
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After an integration by parts, C follows, 



1 

2m 

1 



lUvfrr,- 1 ^ 



1 



^(iha- l d v l 2x2 --gW"<r a )) 



{-iha- 1 !! v l 2x2 - \gWZ<r a ))V 



^f&{Wl* h ){WZ* e )% (56) 

where the arrows above the partials indicate in which direction the derivative is taken 
and the non-Abelian gauge field 

-gW^CTa = maap, (57) 

was introduced with g = h so that W°, has dimensions of ma/h. The SU(2) gauge 
field is then proportional to the SO strength. The contraction on the internal index 
a = x,y, z is understood. From this Lagrangian, the spin current follows from 
derivatives with respect to the gauge field components, 

K = ^ (58) 

which has the correct dimension of a spin current. A first contribution Sm to the 
spin current follows from the second term at the r.h.s. of equation (1561) . 

Jml = - 7^^(-T; h < J a)(-ifia- 1 d v l2x2 - maa p )^ 
v "p 2m 2 

- — ^(iha- 1 d v l 2 x2 - maer p ){--ha a )V, (59) 

2m 2 

and can be explicitly calculated for all three spin components, 

.7(1)" = ~ -^-[tfV*^* - (d v &)a a V] - ha& costp*, (60) 

v "c ima 2 

J(i)l = - j^[&* v d v * - (d v tf)tr v V] - -/k^sin^, (61) 

^ = " ^ [ * tff '^* _ ^*V.*]- (62) 

The "diacolor" contribution (linear in a) mentioned above appears only for the a; and 
y spin polarizations. There is nevertheless another term in the Lagrangian, the third 
term of the r.h.s. in equation (j56[) which plays the role of a gauge symmetry breaking 
term (or a mass term). It brings another contribution J< 2 \ to the spin current, 

7 (2) ; = ^3 2 * f K«>6 + W»oia„) *, (63) 

= ~ha&{a a ,a p }$, (64) 

so 

J i2) ^ = ha^cos^, (65) 

J (2) £ = ~7k*tft sin ^ ( 66 ) 

*7(< = 0, (67) 

that exactly cancels the diacolor contribution in such a way that for all three spin 
polarizations only the "paracolor" contributions survive. Eventually the total spin 
current density, given by 

j; = -^[&<r a d v * - (d v &)<ran (68) 
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is in agreement with the explicit calculations presented above for the z— polarization. 
In the general case, due to the presence GSB term, it does not coincide with the usual 
definition of the current in terms of the symmetrized spin- velocity product (I41j) . This 
property is not very well known and deserves special attention. 

5. Path integral on the ring and voltage quantization 

Many questions in the context of mesoscopic rings concern phase effects. When a 
perpendicular magnetic field is enclosed by the ring there occur peculiar conditions to 
ensure a single- valued wave function [42 . The well known Aharonov-Bohm effect is a 
striking consequence of the peculiar phase relations on a multiply connected structure 
such as a ring. In the presence of SO interactions, phase relations are upgraded to 
spinor interferences since the wavefunctions are now two component objects. Special 
attention must be paid when a SO interaction is present, since a new physical effect 
appears i.e. the spin precession. Spin precession is a two component phase evolution 
in contrast to the one component phase evolution in electromagnetism or U(l) gauge 
theory. 

There is a simple qualitative picture, which consists in the interpretation of the 
interaction term as similar to an effective Zeeman term, ~ — 9 ■ B c ff.. In the case of 
the Rashba SO interaction, the effective magnetic field is 

Vfashba = ap v e p (69) 

and the axis around which the spin precesses rotates when the electron moves along 
the ring. 

This condition appears when we consider the SU(2)— phase accumulation by the 
spinor under transport along the loop |43j : 

h- 1 j ^gW a cr a df = ^ J * a dtp a p . (70) 

This expression follows from the gauge formulation of the problem (4151 FH] , and it is 
apparent for example in the Rashba case in equation (f25j) where the spin of the electron 
appears to be minimally coupled to the electric field through the SU( 2)— kinematic 
momentum 

(-iha^ 1 d v )l 2 x2 - mao-p, (71) 

where macr p = hgW a (T a is a non-Abelian SU(2)— gauge field introduced in 

equation fj57j> . When the operator exp ( iftr 1 § \gW a cr a dr J constructed from 

equation (fT0"|) is applied on * = V*tl+)p + VmJ - )p> the condition 

COS^=l 2x2 * + J(T p sill^=2$ = * (72) 

follows to secure single-valued spinors. We thus obtain a quantization condition 

2irmaa 

= 2n x integer. (73) 

n 

In the case of the standard SO interaction, a = \e\HE/2m 2 c 2 . This condition may 
be rewritten in terms of a typical voltage in the problem defined by V = 2iraE (take 
care here, the electrons on the ring move in a constant potential and the gate voltage 
applied externally to produce the electric field is a different quantity V ga t e ), 

V 

— = integer, (74) 

47rmc z /e 

where one can introduce a quantum of voltage, Vq — Aitmc 2 /e. 
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The same result also follows from more general arguments based on the path 
integral formulation on the ring. As the student has learned in graduate quantum 
mechanics, path integral allow for an alternative and illuminating derivation of 
quantum mechanics generalizing the action principle of classical mechanics. We will 
close our description of ring physics in the presence of SO interactions by deriving in 
detail its path integral and making contact with the quantization conditions derived 
above. 

Let us consider again the Hamiltonian (|25|) and define the eigenstates \n) of the 
operator id v , id v \n) = — n\n), such that periodicity along the ring and normalization 
are satisfied, (ip + 2x\ri) — (tp\n) and L dip\(ip\n)\ 2 — 1. We have 

(<p\n) = (2 7 r)- 1 / 2 e mv , n e Z. (75) 

In order to construct the amplitude for an electron to go from a particular state on 
the ring \ip, a; t) to another similar state at different time \ip' , a'; t'}, we construct the 
path integral in discretized time as 

JV-1 

A(<p,<r,t -+<fJ,o , \t) = Y[(^ +l7 a t+1 \e-' AtIi / h \ Vl ,a t ), (76) 

2=1 

with At = £j_|_i — ti- In the limit At —¥ 0, the evolution operator is expanded to linear 
order in At and one has to evaluate matrix elements of the Hamiltonian between spin 
states. Expanding on the "plane wave" basis, one has 

id v U X 2\<P,<r) = ^(-n)(2 7 r)- 1 / 2 e-^l 2x2 |n,a), (77) 

n,GZ 

and the matrix element of the Hamiltonian reads as 

(pVIHI^) = ^ E h nH °'°' ~ 2^KW]e m( ^ } . (78) 

dlllbUb Zlt\ lb 

ngZ 

When we pass to the evolution operator, one has to evaluate 

7At 7At 

(¥>V|l 2x2 - ■=-B\<p,<r) = WW„, a -<^>'|I%,a), (79) 

n n 

where we make use of the property (<p'\<p) = Xmez 2^r e " ~ m 0T( ^ eT to factorize 
out the term e m ' v ~ v \ and we exponentiate again the Hamiltonian matrix element to 
get 



iAt , , v-^ 1 



nGZ 



(<ff,0 , \l 2 x2-=rH\<p,0) = > 7^ e~^? At[n l2x2 " 2n — ^ e W<P-V). ( 80 ) 



At[n 2 l 2X 2-2n I 



a, a' 



Completing the square and introducing the classical variable v v = a(ip' — if)/ At, one 
eventually arrives at 



-l^m(v v ,l 2X 2-a<'p) 2 ]At 



(<p t+u a t+1 \e- t * ttl " l \<p x ,<T l )= const e 7tl 3 "H<V-^-^ ^ ) (gl) 

L J a,a f 

such that extended to the whole path, we get the symbolic expression 

A(<p,a;t-^(f',a';t')= [ Vtp(t)T [e* f ^(^U^-a^fdt] ^ ^ 

J L Jct.o-' 

where T is a super-operator ordering chronologically all operator products. We note 
that the path integral representation is only performed at the level of the space 
degrees of freedom and not for the spin variables for which we still have an evolution 
operator. Introducing spin coherent states, one could formulate path integrals where 
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L = ^m(v v l 2X 2 -acr p ) 2 . (83) 



both coordinate and spin would be classical variables 45 . From the path integral, we 
can read a "classical Lagrangian" 

1 

2' 

When we expand the square, we get three contributions to the spin precession. If we 
consider the "phase accumulation" for a classical trajectory along a closed path over 
the ring (one round trip), the first term is simply the dynamical phase, which is path 
dependent, 

T / -mv v adipd a , al = ^—0~„>. (84) 

hj 2 n 

It corresponds to the analogue of the optical phase. The second term, the non-Abelian 
equivalent to the Aharonov-Bohm phase, is the path independent (for closed paths) 
contribution discussed in the beginning of this section, 

1-2-K 



h Jo 

or, using a and V defined above, 
V 



1 f *■ , , ,, 2itmu(cT p ) a ^a 

- mv ip a{(Tp) l j y< j>dt= , (85) 

d above, 

2^^^'- (86) 

This is the counterpart of the famous $/<I>o phase in U(l) gauge theory. The third 
term is an additional path-dependent contribution to the phase which takes its origin 
in the initial GSB term of equation ([56]) , 

1 f 27T 1 2 , 2n .. irma 2 5 a ,cT>a 

U 2 ma(g ^' dt= fa" ■ (87) 

This last term produces a phase variation which varies quadratically with the gate 
voltage applied. 

Let us note the particular role of the GSB term in the Hamiltonian. In equation 
(|82p. the transition amplitude is given by a path integral where the action involves a 
term quadratic in the gauge field. In the simpler case of a particle submitted to the 
/7(l)-electromagnetic gauge field, the corresponding transition amplitude would read 
as 

A(<p; t -> ip'; t') = / Vip{t)ei J{h™l+™* A v) d \ 



i.e. it would not include any quadratic term in the U(l) gauge held. The presence of 
such a quadratic term and its contribution (J87I) to the phase accumulation is a direct 
consequence of the GSB term in the SO case. 

6. Summary and Conclusions 

We have presented an overview of the analytical treatment and role of spin-orbit 
interactions in the problem of mesoscopic ideal rings. For the benefit of the student, we 
have placed particular importance to the subtleties of deriving the correct Hamiltonian 
to the problem when posing it in terms of cylindrical coordinates. Such details 
had been overlooked in the literature for some time until reference [33] clarified this 
point. Besides the derivation given in the latter reference we have presented two other 
appealing approaches that may offer a simpler procedure for more involved problems. 
Once the correct Hamiltonian was posed we set out to explicitly derive the ground 
state properties in terms of the charge and spin currents as a function of the spin-orbit 
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strength, pointing out the symmetries of the problem, and making such symmetries 
explicit both for the wavefunctions and the spectrum. Careful attention was paid to 
the correct limiting behavior as the spin-orbit interaction was sent to zero, the level 
degeneracies and all polarization components of the spin current for the lowest filling 
of the ring. 

Using the Lagrangian formulation of the problems we formally derived the 
equilibrium currents, making a connection with recently reported color currents |40) . 
Such an approach also permitted assessing the role of the naturally ocurring gauge 
symmetry breaking term in the Rashba Hamiltonian, that is absent from the intuitive 
formulation of the current in terms of the velocity-spin anticommutator. The results 
point to a physical distinction of the gauge symmetry breaking contribution in the 
X, y polarization components of the spin current. 

We finally addressed topological considerations on the spin-orbit ring. A voltage 
quantization condition was derived due to uniqueness of the wavefunction around the 
ring. A path integral approach to the problem was also formulated in order to formally 
derive a quadratic in spin-orbit strength contribution to the geometrical phase around 
the ring. 

Every one of the topics discussed has been an opportunity to convey to the 
student the subtleties of the quantum formulation of the problem of spin transport on 
mesoscopic rings, detailing the calculations involved and making emphasis on physical 
insight, and not just bare technicalities. 

There are many roads that begin from the material treated here. It is interesting 
to point out that the exact solutions derived are valid both for the limits: i) Slow 
rotation of the Rashba magnetic field with respect to a rapid precession of the spin 
and ii) Rapid rotation of the Rashba field in relation to the precession of the spin. The 
first is the adiabatic limit where Berry phase effects are dominant, while the latter 
relates to the sudden perturbation limit. Understanding physically the full range of 
such behaviors in the presence of a confining potential introduces the physics of lateral 
subbands that may render many interesting effects for mesoscopic spin-orbit rings as 
quantum circuit elements. 
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